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Abstract 



We consider the impurity-induced polariton local states in a dipole-active 
medium. These states present the local optical vibrations which are coherently 
coupled with the induced electromagnetic field. We show that the crossover 
between the polariton and phonon local states takes place within the relativis- 
tically narrow interval near the bottom of the polariton gap. However, the 
resonant phonon-photon interaction leads to the absence of the lower thresh- 
old of the impurity strength. These results are due to the singularity of the 
density of the polariton states, which is generic for any isotropic dipole-active 
phonon mode with the negative dispersion. We show that a weak anisotropy 
removes this singularity and sets a finite threshold for the local polariton 
states. 
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I. INTRODUCTION 



The concept of the polariton local states was for the first time introduced in Refs. . 
Considering a two-level impurity atom interacting with electromagnetic field in an isotropic 
frequency-dispersive medium, the authors discovered the photon-atom bound states. These 
states describe an excited atom which is coherently coupled with the medium polaritons via 
the dipole interaction with the field. If the intra-atom transition frequency is located inside 
the polariton gap, the radiative relaxation of the bound states is suppressed, and the induced 
electromagnetic field is localized around the impurity. It was shown that the two-level model 
can be exactly solved by means of the Bethe-ansatz technic. Authors constructed a complete 
set of the scattering and bound states and analyzed a formation of the polariton-impurity 
band in a system of spatially correlated impurities. 

The results obtained in the above mentioned papers give a basis for the theory of the 
photon localization in the crystals doped with the dipole-active impurities. However, some 
questions need to be clarified. The optical activity of the considered impurity atom was 
associated with its electronic structure. In the case of the "large atomic orbitals", when 
the characteristic scales of the polariton-atom bound states (coherence length, localization 
radius) are comparable with the orbital size, the non-local electronic structure of a defect 
should affect the results. In the opposite case of the "small orbitals", the effects of the 
microscopic anisotropy of the medium need to be taken into account. 

Moreover, as it was shown in Ref. the spatial dispersion of the medium has a sur- 
prisingly drastic effect on the polariton local states. Considering a dipole-active medium 
with a point-like atomic defect (an "electronically rigid" impurity ion, an isotope impurity, 
etc.), authors found a new type of the defect-induced polariton local states. These states 
are composed from phonons and photons localized around the defect, and they are related 
to the phonon local states However, unlike for pure phonon states, there is no lower 

critical value of the impurity strength that needs to be exceeded for the polariton local state 
to arise. It was shown that near the bottom of the polariton gap the localization radius is 
macroscopically large, what supports the long wavelength approximation used in 0. 

In Ref. [0 , we analyzed the polariton local states associated with non-isotope impurities. 
Considering a cubic crystal with a defect which locally affects the crystal elastic constants, 
we found several local states of different parity. All new states arise at the bottom of the 
polariton gap upon infinitesimally small variations of the impurity parameters. This feature 
of the polariton local states was attributed to the singularity of the density of states at the 
bottom of the gap. It was shown that this singularity is caused by the long-wavelength 
polaritons, what provides applicability of the long-wavelength approximation to the local 
states near the gap bottom. 

In the present paper we investigate the crossover between the polariton and phonon 
local states. The long-wavelength nature of the polariton states allows us to analyze the 
crossover within the continuum model of a crystal. In this approximation all details of 
ion displacements within elementary cells are ignored and only resultant dipole momenta 
of the cells are considered. It leaves the polarization and the electromagnetic fields as the 
only relevant dynamical variables. Such approximation is appropriate when the localization 
radius of the considered states greatly exceeds the lattice parameter. 

Our analysis shows that, in agreement with Ref. 0, a single impurity embedded in 
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the medium give rise to the local states in the polariton gap. However, the frequency 
range available for the local states is narrowed by the longitudinal band. We show that 
the states near the bottom of the polariton gap are composed from the long-wavelength 
polaritons. The typical momentum of these polaritons, /Cmax, defines the coherence length 
of the local states, Zcoh ~ aP~^^^. The separation of the local states from the gap bottom 
defines their localization radius, l\oc ~ af3~^, where a is the lattice parameter, and j3 = v/c 
is the ratio between the phonon's and photon's speeds. Despite the atomic size of a defect, 
both characteristic lengths are macroscopic. 

The long-wavelength nature of the polariton states allows us to analyze the crossover 
between the polariton and phonon local states within the continuum approximation. We 
show that the crossover takes place in the relativistically narrow interval, Ai ~ P^q, near 
the bottom of the polariton gap. However, despite the small size of the crossover region, 
the polariton resonance affects the selection of the defects allowing local states in the entire 
gap and leads to the absence of the lower threshold for the impurity strength. This result 
is caused by the singularity of the density of the polariton states, which is generic for any 
dipole-active phonon mode with negative dispersion and isotropic spectrum. We showed 
that a weak crystal anisotropy removes the singularity from the bottom and sets a finite 
threshold for the local polariton states. 



II. POLARITON LOCAL STATES 

Let us consider a dipole-active medium with an embedded impurity. Dynamical equations 
of the medium can be introduced phenomenologically or derived from the microscopic lattice 
equations. In the latter case the displacements of ions within each elementary cell must be 
first expressed via the dipole momentum of a cell, the displacement of its center of mass and 
other similar variables, then the continuum limit of the equations must be worked out. The 
resultant theory describes the long-wavelength dynamics of a crystal where the electric field 
effectively interacts with the optical phonons only: 

(c.^-i7^)a = -^E. + a^no), (1) 



{u^ - c^k'^) Ek = -47r (a;^ - c^k"^) P\\ + uj'^P ^ . (2) 

Here 7\ and Ek are Fourier components of the polarization and electric fields, V{0) is 
the defect polarization, (k) is the dynamical matrix of the medium, Py and P^ are 
the longitudinal and the transverse projectors in the momentum space, a is the lattice 
parameter, V is the sample volume, and (P is the photon-phonon coupling parameter. 

Eq.(l) imphes that the defect is electrically identical with the replaced host ion. In this 
case the "defect strength" can be expressed as a = u"^5^/^ — ijj'^5jji/jji, where ji and 7 are the 
reduced mass of the elementary cell and the elastic constant of the nearest-neighbor bonds, 

and ^7 are their impurity-induced variations, and a; ^ is the characteristic frequency 
depending on the properties of a defect in the crystal. 

In the isotropic medium the dynamical matrix can be presented as follows: 
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= QfP|l +QiPx, (3) 

where r2| and are the frequencies of the longitudinal and transverse phonons, respectively. 
Below we assume a negative dispersion in the phonon branches; so they have the following 
long- wavelength asymptote: 

Ql (k) vl k\ (4) 

nlik) ^nl-v^e, (4') 

where the parameters vj_ and v\\ set the ranges of the typical phonon velocities, and Qq is 
the phonon activation frequency. For the "order of the magnitude" estimates we assume 
v± ~ v\\ ~ lO^m/s, Qq ~ v/a ~ d, and a ~ a Ql, where a is a numerical parameter. 

Solving Eqs.(l-3), one can express Ek and Vk in terms of the defect polarization, 7^(0), 
and obtain the spectral equation: 



3 \27tJ 



dk 



ijp- — c^k^ 



aliu'), (5) 



where the integration is extended over the first Brillouin zone. This equation defines spectra 
of all, extended and local, excitations. The extended states form continuous bands and their 
dispersion relations are determined by the poles of the integrand in Eq.(5). In the isotropic 
medium there is a single longitudinal branch, o;^ = Q| + d^, whereas the transverse band 
contains two polariton branches [Fig.l]: 



Q± (k) = \ 



\l (Qx + Ckf + Ci2 ± ^J{n^-Ckf + 



(6) 



The lower branch, {k) , is activationless and non-monotonic. Analysis of its "large- 
momenta" asymptote, 

Q_ {k) ^ (k) (i - ^) , (7) 

shows that the lower branch reaches its maximum at the point 

1/2 



(8) 



Since f/Cmax ~ (^^^"^^0 -C VLq and c/cmax ~ /3~^^^f2o ^ the maximum of Vl?_{k) is 
located in the long-wavelength region but far away from the cross-resonance point. Using 
Eqs.(7-8) one can find: 

^Ql-i^^^nl- 2pdQo. (9) 

The polariton branches are separated by the gap which extends from l^^^x tof^Q + d^. 
However, because the longitudinal branch overlaps the top part of the polariton gap, the 
true spectral gap, with no modes of any kind inside, is between O^^^ and the minimum of 
the longitudinal branch, 0,^^^^ — fll + d"^ — Ay. The spectral gap is open only if the width of 
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the polariton gap, (P + Ai, exceeds the width of the longitudinal band Ay. Further below 
we assume this and consider only local states inside the gap. 

The frequency of the local state can be found from Eq.(5). The function /(cu^) in this 
equation presents a sum of the "longitudinal" and "transverse" terms, and I±{uj'^). 

The transverse integral can be written in the form: 



When the frequency approaches flmax, this integral diverges at the surface = k^^^. Near 
the bottom of the gap {uj > l^max) j the transverse term dominates in /(u;^) and the region 
of /c ~ /Cmax givcs the major contribution to /^(cu^). As it follows from Eqs.(4, 9), in this 
region we can approximate: 

^-{k)^^ls.-'^vl {k-k^^f. (11) 

It allows us to evaluate /(o;^) in Eq.(5) and obtain the frequency of the local state: 

V^' - ^leo. ~ ^ {ak^s^f - a(3no. (12) 

In order to find the localization radius of the state we need to consider the spatial 
distribution of the electric and the polarization fields. Using the solution of Eqs.(l-2), we 
obtain: 



q/q3 I" 

E (r) = -— / dkexp (ikr) 



+ 



UJ 



7^(0). (13) 



Near the bottom of the gap the longitudinal field is small comparing to the transverse term. 
Taking into account that the region of A; ~ /cmax gives the main contribution to E^ (r) and 
employing Eqs. (11-12), one can find: 

(r) ^ 4. ^ {nx [n x P(0)]} "° ''''-^ ^-"''> . (14) 



In a similar way one can obtain: 



^ r \ r r -o/nMi sin (r/Cmax) exp (-rx) 

'P±{r) fti {nx [n x P(0)]} ^ f_ ^ (15) 



>'max 



where n is the unit radial vector, and 



- Ql 



max 



(a/^max)^ ~ otiia ^ (16) 



is the inverse localization radius of the considered state. 

Equations (14-15) allow us to estimate the energy partition between the phonons and 
photons bound in the local state. Taking into account that the polarization and the electric 
fields are concentrated within the x~^-range from the defect, we can evaluate their energies: 
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W^^^n - J V\r)dr ~ V\0)k;,l><-\ (17) 
Wf,,i, ~ / E^(r)dr ~ V'iO) (17') 
It gives us the "order of the magnitude" estimate of the ratio of the energies: 



W field / ^0 \^ o2 



P\ (18) 



III. CROSSOVER BETWEEN THE POLARITON AND PHONON LOCAL 

STATES 

The features of the polariton states near the gap bottom are caused by the singularity 
of the density of states in the lower polariton band. Since fl'^_{k) reaches its maximum at 
the surface of a finite area inside the Brillouin zone, the density of states diverges at ilmax^ 



dnl{k) 
dk 



aiak 



\2 



max ; 



(19) 



This singularity is provided by the long-wavelength {k ~ Knax) polaritons, which also dom- 
inate in the local states near the bottom of the gap. It explains the macroscopic size of the 
coherence length, /coh = ^ a,(3~^^^, and the localization radius, /loc = ~ a(3~^. Also, 
because the typical momentum of these modes is far away from the cross-resonance point 
(^max ko), the phonon content of the local states greatly exceeds their photon content. 
The local polariton states appear first at the bottom of the gap for a = +0 and move inside 
the gap upon increase of the defect strength. This further weakens the photon content of 
the polariton local states and transforms them into the ordinary phonon states. 

To investigate the crossover between the polariton and phonon local states we need to 
analyze Eqs.(l-2,5) in the entire spectral gap. Away from the gap bottom both terms of 
Eq.(5), I\\{uj^) and I±{u^), become comparable and calculations here require knowledge of 
the phonon dispersion laws f^|(/i;) and Q^lk) in the entire Brillouin zone. Such detailed 
information is not consistent with the approximations employed in our model. In the short- 
wavelength region a crystal anisotropy cannot be neglected, and ion displacements within 
elementary cells must be considered in all details. However, assuming that the crossover 
takes place near the gap bottom, we proceed with calculations noting that in the isotropic 
model it is enough to know the densities of states in the phonon bands. We use a simple 
approximation accounting only for Kohn's singularities at the band boundaries: 



Pphon(£) = V - ^min) «ax " (20) 

phon 



where Aphon is the width of the phonon band; a;^;^ and (^^ax band boundaries. 

To calculate /y (o;^) we transform it into the integral over the longitudinal band: 
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II min 

where f^mm = + — Ay. Using Eq.(20) for p\\{£) one can obtain: 
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^11 K) = - ^ l^V^-^- + - - V ^-in - '^^J • (22) 

To evaluate the transverse integral wc present it as a sum of two terms by separating the 
Brillouin zone into two parts: 

The first term, /^(c<j^), involves integration over k < k'^ including the cross-resonance region, 
and presents the integral over the remaining part of the Brillouin zone, k > k'^. The 

separating momentum /cq is convenient to fix by the condition ^^^.(^o) — ^max [Fig-1]- 

Since k'^ is located far away from the cross-resonance point, the product ck in iKou'^) 
greatly exceeds the typical phonon frequencies and we can set there ck — oo: 

Using Eq.(20) for p±{e), we can calculate the last integral: 

1 ^ r 

I'li'^^) = inr- (25 - 1) arcsin ^ - (1 - 5i) + |1 - 6\ F{u;^) , (25) 



37rA 

where 



= \ V^O=^VMp) ~ ~ ° (26) 

-2arctany^|Q for u;^>nl, 

and we denote 6i = {A± — Ai) /A±, and 5 = (A_l + lo'^ — fig) / A^. 

Considering cu^ far away from the gap bottom (a;^ — fl^^ ^ Ai) , one can obtain: 



^uj^ + A^-Ql-^uj^-Ql^ . (27) 



3A2 

In the opposite limit, cu — > flma^, Eqs. (25-26) lead to the logarithmic singularity: 

'"<^^* " S'""^- '''' 

The second term in Eq.(23) includes integration over the cross-resonance region where 
the polariton effects cannot be neglected. Recalling that /i(<^^) presents the integral over 
the long- wavelength modes and using Eqs. (4, 20), we can rewrite I±{uj'^) in the form: 



7 



Ml 

7n2 „ [e - 



{nl - e) p±{e) de 



{e - [e - Qg) + {ujdv/cf SttA^ {x - a+) {x - a_) 



(29) 



Here h — / Ai is a large parameter and two poles of the integrand are given by the 
equation: 



2Ai 



(30) 



where a;^ = JIq ± Ai, so that a;^ coincides with the gap's bottom, l^^ax- 

Direct calculations show that in the interval a;^ < ^ '^+) where a+ and a_ are 
complex- valued, /K<^^) has the form: 



3^/2 - 



3r/V 1 
arctan — 



In 



2r] 



V 1 + u 
h arctan 



+ 



(31) 



4z/ (l + z/)2 + r/2_ 
where the parameters u and 77 are defined by the equations: 



2 2 
— ?7 



2Ai 



2 1/77 



a;' 



CO — CO 



2\V2 



2Ai 



(32) 
(32') 



In agreement with the results obtained in the previous section, Eqs. (31-32') give the "square- 
root singularity" of I±{uj'^) at the gap bottom: 



/', (c^^) cx 



8A1 



3A 



3/2 



-1/2 



(33) 



For > f^o + Ai, when a± have real negative values, we obtain: 



1/2 
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3/2 
± 



1 - 



\a+ 


arctan 


a+l 


-1/2 




a_ 


arctan 


a_ 


-1/2 








a_ 





(34) 



Eqs. (22, 25-26, 30-34) allow us to analyze the local states in the entire gap. Since any 
local state presents a superposition of different modes, we can qualitatively describe its 
composition comparing contributions to /(cu^) from different bands. 

As we showed, the function /(cj^) has the singularity at the gap bottom caused by the 
long-wavelength {k ~ fcmax) polaritons. They also give the major contribution to I{uj^) in 
the interval, f^max < ci;^ -C I^q. However, even at a;^ = fl^ their contribution to /(a;^) is 
substantially weakened: 



3Ai 




(35) 
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where the leading term is provided by the short-wavelength {k > k'^) polaritons, which are 
physically indistinguishable from the transverse-optical phonons. 

Outside the interval f^^ax < < ^max + 2Ai the function /(cu^) has basically the 
phonon structure. All terms related to the polariton singularity are weakened there by the 
power factors of the small parameter Ai/A^. At the upper boundary of the gap we obtain: 



I{VL 



2. ) = ^ 



X 2 




(36) 



where x = A± / Ay and y = (P / Ay — 1. 

-^(^min) determines the upper limit of the impurity parameter which allows the local 
state to exist. Its value can be positive or negative depending on the relations between the 
ion plasma frequency d? and the widths of phonon bands A_|_ and Ay. Analysis shows that 
-^(^min) is positive only if the width of the spectral gap does not exceed a certain value, 
- A|| < (A|| / 8) (2 - Ax/A||) . Therefore, only in the ' 'narrow" gaps all local states are 
associated with impurities of the same type. In the case of the isotope defects, for instance, 
the local states arise only in the presence of light impurities. However, if rf^ — Ay exceeds 
the critical value, the local states near the top and the bottom of the gap are associated with 
the different types of defects. For isotopes those are heavy and light impurities, respectively. 
The frequency regions of the corresponding states are limited inside the gap. The upper 
frequency for the "light" states and lower frequency for the "heavy" ones are defined by the 
equations, uj'^I{uj'^) = 1 and /(c^^^) = 0, respectively. 



IV. ANISOTROPY-INDUCED THRESHOLD FOR THE LOCAL STATES 

The considered above impurity-induced local states appear right at the bottom of the 
polariton gap upon infinitesimally small variations of the impurity parameters. This is in 
contrast with 3-D phonon systems where a lower threshold for the local states always exists. 
The general theorem regarding the existence of the threshold for local states in bandgaps of 
the periodic systems was given in Ref. 0]. However, the proof suggests that the density of 
states is regular everywhere. In our model, the absence of the threshold is caused by the sin- 
gularity of the density of states in the lower polariton band, p (o;^) oc (f^max " ^'^) ^^'^ ■ This 
singularity is provided by the long-wavelength polaritons and is generic for any dipole-active 
phonon mode with negative dispersion and isotropic spectrum. 

However, even a weak crystal anisotropy can remove the singularity of the density of 
states from the gap bottom. A dipole-active isotropic medium presents the "zero-order" 
approximation of a cubic crystal. The anisotropic terms of the phonon dispersion law in the 
cubic crystals appear beyond the quadratic approximation: 

ni(k) = nl- vie + x^l {akf F (^k) (37) 

where x is the small parameter, and F ^kj is some anisotropic function. The anisotropic 

term, as one can see from Eq.(7), makes the position ( fcmax) and the value of the maxima 
in the lower polariton branch ( f^^ax) dependent on the crystallographic direction. In the 
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case of a weak anisotropy, the first effect can be neglected, what leads to the following 
approximation of fi^ (k) near the surface k^ = 



.2 

''max ■ 



(k) ^ ^]Lx [l + XF (kj (aA:^ax)'J - {k - k^^^f 
The small magnitude of the anisotropic term, 



2 

max' 



XF [^kj (a/^max)"^ 

allows us to evaluate the asymptote of the density of states: 



(38) 



(39) 



(40) 



Since the crystal anisotropy destroys the continuous degeneration of the polariton spectrum, 
it removes the singularity of the density of states from the gap bottom. The anisotropy- 
induced shift of the singularity. 



F k 



(41) 



where 



F k 



1 is the "average" value of the anisotropic function, defines a finite 



threshold for the polariton local states: 



(42) 



V. SUMMARY AND CONCLUSIONS 

We considered the local states in an isotropic dipole-active medium. It was shown that, 
in agreement with Ref. a single impurity embedded in the medium give rise to a local 
state in the polariton gap. The frequency interval available for the local states is extended 
from the top of the lower polariton band, f^^ax' ^'^ bottom of the longitudinal band. 

The polariton local states appear right at the bottom of the polariton gap for an in- 
finitesimally small value of the impurity strength. Our analysis showed that the local states 
near f^^^^^ are composed from the long-wavelength polaritons. The typical momentum of 
these polaritons, kynux, defines the coherence length of the local states, /coh ~ o (3'^/"^. The 
separation of the local states from the gap bottom, ^Juj"^ — fi^^x; defines their localization 
radius, l\oc ~ a(3~^a~^. Despite the atomic size of a defect, both characteristic lengths are 
macroscopic near the gap bottom. 

These results, and the absence of the lower threshold for the impurity strength, are due 
to the singularity of the density of the polariton states, p_ (a;^) oc (fi^ax ~ ^'^) • This 
singularity is caused by the long-wavelength polaritons and is generic for any dipole-active 
phonon mode with negative dispersion and isotropic spectrum. We showed that a weak 
crystal anisotropy removes the singularity from the bottom, p_ (cj^) oc (^^^ax — + ^) , 
and sets a finite threshold for the local polariton states, amin ~ A^/"^. 
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The long-wavelength nature of the polariton states allowed us to analyze the crossover 
between the polariton and phonon local states within the continuum approximation. We 
found that the crossover takes place in the relativistically narrow interval, Ai ~ f^^o, near 
the bottom of the polariton gap. However, despite the small size of the crossover region, 
the polariton resonance affects the selection of the defects allowing local states in the entire 
gap. We showed that in the "narrow" gaps there must be only light isotope impurities. 

In the presence of a macroscopic concentration of impurities, the local states can provide 
the transmission inside the spectral gap. Since l\oc greatly exceeds /coh, the transmission 
regime critically depends on the impurity concentration, n. When n~^^^ <^ /loc the proba- 
bility of the tunneling of excitations from one impurity to another is exponentially small, 
what corresponds to the localized regime. For /coh < n~^/^ < hoc the resonant tunelling 
gives rise to the diffuse propagation of the radiation. When n~^^^ becomes less than /coh, 
the polariton-impurity band begins to form, and the transmission regime regains properties 
of the coherent propagation. 

The local polariton states provide a "phonon assisted locahzation" of electromagnetic 
waves. However, our estimate of the energy partition shows a strong suppression of the 
photon content of these states, Wfieid/Wmech ~ P^- It is caused by the fact that /cmax is 
much greater than the cross-resonance momentum /cq. To eliminate this disproportion, we 
need to lower the group velocity of electromagnetic waves in the active medium. It can 
be achieved in the medium placed inside a waveguide. For instance, the dispersion law 
of the modes propagating in the parallel- plate waveguide, oj/c = provides 
the reduction of their group velocity in the long-wavelength region. Estimates show that 
for / ~ lO^a the phonon and photon velocities become comparable in the cross-resonance 
region. The detailed analysis of the possibility to obtain the local modes in a thin waveguide 
filled with a dipole-active medium will be presented elsewhere. 
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